% model = createpde;
% geometryFromEdges(model,@lshapeg);
% mesh = generateMesh(model);
% pdeplot(model)
% a = mesh.Nodes;
% b = mesh.Elements;
% 
% load('test.mat')
% nodes = [x;y];
% elements = double(tri)'+1;
% figure
% pdeplot(nodes,elements,"XYData",qh(:,1),"XYStyle","flat") 
% figure
% pdeplot(nodes,elements,"XYData",qh(:,1),"ZData",qh(:,1),"XYStyle","flat", ...
%     "ZStyle","discontinuous") 

load('test2.mat')
nodes = [x;y];
elements = double(tri)'+1;
figure
pdeplot(nodes,elements,"XYData",qh(:,1),"XYStyle","flat") 
figure
pdeplot(nodes,elements,"XYData",qh(:,1),"ZData",qh(:,1),"XYStyle","flat", ...
    "ZStyle","discontinuous") 